Draft version July 20, 2010 

Preprint typeset using L£-T^X style cmulatcapj v. 6/22/04 



O 

00 



Oh- 
i < 

o 

o3 



> 
oo 
in 
o 

O 
O 



13 



A BIPOLAR OUTFLOW FROM THE MASSIVE PROTOSTELLAR CORE W51E2-E 

Hui Shi 1 , Jun-Hui Zhao 2 , J. L. Han 1 

Draft version July 20, 2010 

ABSTRACT 

We present high resolution images of the bipolar outflow from W51e2, which are produced from the 
Submillimeter Array archival data observed for CO (3-2) and HCN(4-3) lines with angular resolutions 
of 0."8 x 0."6 and 0."3 x 0."2, respectively. The images show that the powerful outflow originates 
from the protostellar core W51e2-E rather than from the ultracompact Hn region W51e2-W. The 
kinematic timescale of the outflow from W51e2-E is about 1000 yr, younger than the age (~5000yr) 
of the ultracompact Hn region W51e2-W. A large mass-loss rate of ~ 1 x 10 _3 Mq yr -1 and a high 
mechanical power of 120 Lq are inferred, suggesting that an O star or a cluster of B stars are forming 
in W51e2-E. The observed outflow activity along with the inferred large accretion rate indicates that 
at present W51e2-E is in a rapid phase of star formation. 

Subject headings: ISM: individual objects (W51e2) — ISM: jets and outflows — stars: formation 



1. INTRODUCTION 

Molecular outflows from protostellar cores provide a 
critical means to transport angular momentum from the 
accretion disk into the surrounding environment. The 
molecular outflows and infalls vigorously affect the tur- 
bulence and dissipation of molecular gas in molecular 
cores, playing an important role in the formation and 
evolution of massive stars. Observations of outflows can 
reveal the history of mass-loss processes in a protostellar 
system. In particular, high angular resolution observa- 
tions of outflows can precisely determine the origin of 
an outflow and separate the outflow from the infall in a 
protostellar core. Nearly half of the observed molecular 
outflows are driven by massive protostars (Lboi > 10 3 £©) 
and have a typical dynamical age of ~ 10 4 yr and a ma ss 
outflow rate of ~ 10 -3 M© yr" 1 (e.g. IChurchwellll2002ft 

W5 1e2, located at a distance of 5.1 kpc (jXu et al.l 
2009), is a prototypical core for massive star formation. 
The molecular core is associated with an ult racompact 
Hn region (e.g. Scott 1978; Gau me~et al.l!993l) and a pos- 



sible inflow (e.g. |Ho fc Youndfl 996: Zh ang et al.l 119981 : 
ISollins et al]|2004 IShi et al.ll2010D . With high-resolution 
observations using the Submillimeter Array (SMA) 3 at 
the wavelengths of .85 and 1.3 mm , W51e2 was resolved 
into four sub-cores (|Shi et al.ll2010l) including the ultra- 
compact Hn region (W51e2-W) and a massive protostel- 
lar core (W51e2 -E). From the a nalysis of the HCN(4-3) 
absorption line, IShi et ail (j2010l ) showed that the bright 
dust core, W51e2-E, 0."9 east of W51e2-W, dominates 
the mass accretion in the region. A bi polar outflow has 
been detected in the W51e2 region (|Keto fc Klaassenl 
120081) based on the SMA observations of CO (2-1) line 
with an angular resolution of ~ 1". However, the reso- 
lution of the CO (2-1) observation was not adequate to 
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determine the origin of the molecular outflow. 

High angular resolution observations are imperative to 
identify the origin of the molecular outflow. In this Let- 
ter, we present high-resolution images of the bipolar out- 
flow using the SMA data of the HCN(4-3) and CO (3-2) 
lines observed at 0.85 and 0.87mm with angular reso- 
lutions of 0."3 x 0."2 and 0."8 x 0."6, respectively. We 
determine and discuss the physical properties of W51e2- 
E and W51e2-W cores, and assess the roles of the outflow 
in the process of massive star formation in W51e2. 

2. DATA REDUCTION 

The interferometer data for the HCN(4-3) and CO (3- 
2) lines at v Q = 354.505 and 345.796 GHz were ac- 
quired from the SMA archive, observed on 2007 June 
18 and 2008 July 13 with the "very extended" and "ex- 
tended" array configurations, respectively. The reduc- 
tion for the HCN (4-3) line data was made in Miriad 
(jSault et al.1 119951 ) following the reduction in structions 
for SM A data 4 , which has been discussed in IShi et all 
(2010). We constructed image cubes for the HCN(4-3) 
line with a channel width of 1 km s _1 and an FWHM 
beam of 0."33 x 0."24 (P.A.=68°). The typical rms noise 
in a channel image is 0.07 Jy beam -1 . 

The CO (3-2) line is included in the SMA observations 
at 0.87 mm f or the polari z ation study of W51e2 that was 
discussed by I Tang et al.l (|2009D . The cross-polarization 
data were flagged before processing the data in Miriad. 
The data scans on 3C454.3 were used for bandpass cal- 
ibration. QSO J1751 + 096 was used for calibration of 
the complex gains, and Uranus was used to determine the 
flux-density scale. The continuum level was subtracted 
from the line channels using UVLIN. The continuum-free 
line data were used for imaging with the natural weight- 
ing and a channel width of 1 km s _1 . A typical rms noise 
of 0.14 Jy beam" 1 and an FWHM beam of 0."79 x 0."62 
(P.A.=— 80°) were achieved for a channel image. 

The shortest UV lengths in the visibility data are 30 
kA and 54 kA, corresponding to angular sizes of about 7" 
and 4" for the CO and HCN structure, respectively. The 
compact structure (<3") of the outflow in both CO and 

4 \http://www. cfa. harvard. edu/sma/miriad\ 
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Fig. 1. — (a) SMA image of W51e2 for CO(3-2) line emission at 0.87mm, integrated from the velocity ranges of —124 to —12 km s -1 
(blue contours) and +10 to +116 km s -1 (red contours). The integrated absorption (dashed contours) corresponds to the velocity range 
— 11 to +9 km s — 1 . The contours are ±5(J X 2 n Jy beam" 1 km s -1 (n = 0, 1, 2, 3, ...), where a =1.2, 1.3, and 0.63 Jy beam" 1 km s -1 
for the blue-, red-shifted emission, and absorption, respectively. The FWHM beam is 0."79 X 0."62 (P.A.=— 80°) shown at top right. 
The top left inset is the CO(3-2) absorption spectrum toward W51e2-E. The bottom right inset is the intensity-weighted velocity image 
derived from the CO(3-2) line data (the color wedge representing the velocity range —64 to +76 km s -1 relative to Vlsr=53.9 km s —1 ), 
excluding the absorption and insignificant emission (< 3cr). Crosses mark the positions of the sub-cores in W51e2 (Shi ct al. 2010). (b) 
SMA image of W51e2 for HCN(4-3) line at 0.85 mm, integrated from the same velocity range as CO(3-2) in (a). The contours are ±5c X 2 n 
Jy beam -1 km s — 1 (n = 0, 1, 2, 3, ...), where a =0.71, 0.67, and 0.31 Jy beam -1 km s -1 for the blue-, red-shifted emission, and absorption, 
respectively. The FWHM beam is 0."33 X 0."24 (P.A.=68°) shown at bottom right. 



HCN appears to be adequately sampled in the SMA ob- 
servations. However, the total line flux density of the out- 
flow may be underestimated due to missing short spac- 
ing. From a visibility model for the outflow source with 
a uv coverage identical to the SMA data, we assess that 
<2% and <8% in the CO and HCN line flux densities 
may be underestimated, respectively. 

3. RESULTS 

3.1. Morphologies of the bipolar outflow 

Figure[IJa) shows the CO(3-2) image of W51e2 for the 
integrated blue-shifted (—124 to —12 km s -1 ) and red- 
shifted (+10 to +116 km s^ 1 ) emission and absorption 
(—11 to +9 km s -1 ) with respect to the s ystemic veloc- 
ity of VLsr = 53.9 km s" 1 (|Shi et al.ll20"Tot) . The outflow 
velocities discussed throughout the Letter are with re- 
spect to this systemic velocity. The low velocity limits 
were determined by the velocities at the boundary of the 
absorption region. The high velocity limits correspond 
to the velocities at which the high-velocity wing profiles 
drop to the la level. The ranges of velocities for the 
outflow wings are also valid for the CO(2-l) and HCN(4- 
3) line profiles. The CO (3-2) line profiles appear to be 
contaminated by three narrow (AVfwhm ~ 6 km s" 1 ) 
anonymous lines at the velocities of -104, -89 and +70 
km s -1 . Prior to integrating the CO (3-2) line from the 
outflow, the line emission in the velocity ranges from -117 
to -79 and from +62 to +76 km s" 1 is blanked. Thus, 
about 20% and 10% of the line fluxes are missing in the 
blue- and red-shifted outflow lobes shown in Figure [l|a), 



respectively. 

The CO (3-2) image shows a stronger blue-shifted lobe 
in the southeast with a position angle of 140° ± 6° 
(from north to east) and a weaker red-shifted lobe in 
northwest (—40° ± 11°), which are in agreement with 
the outflow structure revealed in the CO (2-1) image 
(|Keto fc Klaassenll2008f ). With Gaussian fitting, we de- 
termined the intrinsic sizes of 2. "3 x 0."6 and l."l x 0."6 
for the blue- and red-shifted lobes, respectively. Our high 
resolution CO (3-2) image clearly shows that the molec- 
ular bipolar outflow originates from W51e2-E (the pro- 
tostellar core) rather than W51e2-W (the ultracompact 
Hh region). The bright emission ridge along the north- 
east edge of the blue-shifted lobe corresponds to a high- 
velocity component of the outflow, or an outflow jet (see 
the velocity distribution image in the lower right inset 
of Figure QJa)). The weak, relatively extended emission 
along the southwest edge of the blue-shifted lobe may 
indicate a backflow of the molecular gas from the inter- 
action between the outflows and the medium surround- 
ing the W51e2 core. In addition, a velocity gradient is 
present along the major axis of the red-shifted lobe, with 
higher velocities near the core, indicating that the gas 
in the red-shifted outflow is decelerated. By averaging 
the intensity-weighted velocities along the bright emis- 
sion ridges in both the blue- and red-shifted lobes, we 
infer the mean radial velocities of Vbiuc = — 34 ± 7 and 
Kcd = +52 ± 2 km s^ 1 for the outflow lobes. The CO(3- 
2) spectrum toward W51e2-E also shows a large fraction 
of absorption to be red-shifted relative to the systemic 
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velocity (top left inset of Figure QJa)), suggesting the 
presence of infalling gas, which agrees with the results 
from the observ ations of the CO(2-l) and HCN(4-3) lines 
(|Shi et al.ll2010h . 

The HCN(4-3) image in Figure W[b) shows the high- 
density gas components in the bipolar outflow. At the 
resolution of 0."3 x 0."2, the SMA image unambiguously 
verifies that W51e2-E is the driving source for the bipo- 
lar outflow in the region. We also noticed that the HCN 
emission traces the northeast edge of the blue-shifted 
lobe where the CO (3-2) emission shows a sharp gradi- 
ent both in intensity and velocity, suggesting that the 
line emission of the molecular HCN is likely enhanced 
near the shocked or compressed r egions of an outflow 
l|Bachiller fc Perez Gutierredll997t) . In addition, at the 
southeastern tip of the blue-shifted outflow lobe, an "L" 
shape structure is detected in the HCN(4-3) emission 
with a double spectral feature at V w —14 and +16 km 
s _1 , each of the spectral profiles having a line width of 
AVfwhm ~ 15 km s _1 . The HCN line emission suggests 
that the gas in the region is swept up by the suspected 
bow shock at the tip of the outflow. 

The polarization emission detected in the W51e2 re- 
gion shows a "hourglass" shape in magnetic fields, which 
is probably relate d to the gas infall and disk rotation 
(|Tang et al.|[20Q9l ). Our images for the bipolar outflow 
in both the CO(3-2) and HCN(4-3) lines (Figures [[]» 
and (b)) show that, located near the "hourglass" center, 
W51e2-E appears to be the energy source responsible for 
both the magnetic fields and the bipolar outflow. The 
"hourglass" structure appears to slightly extend along 
the major axis of the outflow, suggesting that the mag- 
netic fields might also be tangled with the outflow mo- 
tions. 

3.2. Physical properties 

The bipolar outflows from W51e2-E have been de- 
tected in both CO(3-2) and CO(2-l) with the SMA. In 
order to assess the physical conditions of the outflow 
gas, we re-processed the CO(2-l) data and convolved 
the CO (3-2) image to the same beam size as that of 
the CO(2-l) image (l."4 x 0."7). For both the CO(3- 
2) and CO(2-l) lines, the line peak (v — v v ) brightness 
temperatures ATl(v p ) of the blue- and red-shifted line 
profiles can be determined from the maximum positions 
of the CO (3-2) line fluxes in Figured] Then, the optical 
depth at the line peak can be es timated if the gas is in 
local thermal equilibrium (LTE) (jRohlfs & Wilsonll2004ft 
and the cosmic microwave background (CMB) radiation 
is negligible, 

AT L (u) / aa 
T (u) \ ° x 

where T cx is the excitation temperature of the transition 
andTo(^) = hu/k. Thus, each of r(^)( 3 _ 2 ) and t(v)(2-i) 
can be determined independently from Equation (fTJ) at a 
given T ex . On the other hand, under the assumption of 
LTE, the optical-depth ratio between CO (3-2) and CO (2- 
1) can be independently expressed as 

^0(3-2) 



T L {v) 



-In 



1 



1 



(1) 



TL(v) (3- 2) _ (S> )(3-2) e" 
T L{v)(2- 



1 



T (3-3) 



OV) 



(2-1) e 



1 



(2) 



where S and fi are the line strength and the permanent 
electric dipole moment, respectively. For the two given 



transitions, such as CO(3-2) and CO(2-l), the ratio of the 
optical depth {tl(s-2) I t l{2-i)) depends only on the exci- 
tation temperature T cx (Equation (|2"]l). By giving an ini- 
tial value of T ex , the three unknown quantities (t£( 3 _ 2 ), 
7i(2-i), and T ex ) can be solved with a few iterative steps 
from Equations ([T]) and {2J , constrained by the measured 
quantities ATW 3 _ 2 ) an °- AT i ( 2 -i) for each of the blue- 
and red-shifted lobes. 

Table [1] summarizes the results of the physical parame- 
ters determined for the molecular outflow. For the quan- 
tities directly determined from the CO images, the esti- 
mated lcr uncertainties are given. The excitation tem- 
peratures of 120 and 65 K are derived for the blue- and 
red-shifted CO outflow lobes, respectively. The optical 
depths of the blue-shifted CO lines are T£( 3 _ 2 ) =3.9 and 
t l(2-i) — 1-9, and the red-shifted lobe are 



(3-2) 



=1.3 



and 7i(2-i) = 0-7- Uncertainties of ^30% in the derived 
values for T ex and tl are mainly due to the uncertainties 
in the determination of the line flux densities. The H 2 
column density (6.4xl0 22 cm~ 2 ) in the blue-shifted lobe 
appears to be a few times greater than the value derived 
from the red-shifted lobe (1.3 x 10 22 cm -2 ). Considering 
the intrinsic size of the outflow, we infer that the total 
outflow mass is about 1.4 M Q . From the intrinsic _sizes 
(#Maj) and the intensity- weighted mean velocity (V) of 
the outflow components, the kinematic timescale (£ ut) 
is found to be ~1000 yr if the outflow inclination angle is 
taken as i ~ 45°. Therefore, the mass- loss rate and mo- 
mentum rate are lxlO -3 M Q yr -1 and 0.04 M km s _1 
yr -1 from W51e2-E, respectively. The mass-loss rate has 
the same order of magnitude as the value of th e accre- 
tion rate (> 1 x 10~ 3 M Q yr-^ IShi et~aT1l2010D derived 
from the HCN(4-3) absorption line toward W51e2-E. The 
momentum rate corresponds to a mechanical power of 
120 Lq, which is at least an order of magnitude greater 
than that of an e arly type B star (e.g. IChurchwelilll999l : 
lArce et aT]|2007l ). suggesting that the protostellar core 
W51e2-E is forming an O type star or a cluster of B type 
stars. 

4. DISCUSSIONS 

The current star formation activities in the W51e2 
complex are shown by the powerful bipolar outflow orig- 
inating from W51e2-E, in contrast to the ultracompac t 
H n region W51e2-W detected at 6 cm (|Scotti 11978ft . 
The central region of W51e2-W (< 0.06"), character- 
ized by a hypercompact Hn component (see the H26a 
line emission region in Figure ^ with T e = 12, 000 K, is 
surrounded by an expanding ionized compo nent with a 
mean temperature of 4900 K (|Shi et al.ll2010h . The mean 
density of 3.0 x 10 5 cm -3 and a Lyman photon rate of 
3.0 x 10 48 s" 1 were also derived. Using a classic model fo r 
the expansion of an Hn region (jGarav & Liza no 1999), 
we estimated the dynamic age of the Hn region as the 
time of the sound wave traveling from the initial ion- 
izat ion front at a ra dius close to the Stromgren radius 
(r«: lStromgrenin~939f ) to the current ionization front at a 
radius of rs : 



t 



7/4 



1 



(3) 



Here ri can be determined from radio continuum obser- 
vations of the Hn region, we adopted the size of < 2.5" 
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TABLE 1 

Outflow Parameters from CO(3-2) and CO(2-l) Lines. 



Quantities 




CO(3-2) 




CO(2-l) 




Blue 


Red 


Blue 


Red 


V range (km s _1 ) 


-124 to -12 


±10 to ±116 


-124 to -12 


±10 to ±116 


V (km s" 1 ) 


-34 ± 7 


±52 ±2 






AV (km s" 1 ) 3 


19±5 


33±8 


14±2 


19±4 


Outflow intrinsic size #Maj (") 


2.3 ±0.3 


1.1 ± 0.4 






Outflow intrinsic size #Mm (") 


0.6 ±0.5 


0.6 ±0.6 






Outflow P. A. (°) 


140 ± 6 


-40 ± 11 






/ AI(v)dv (Jy beam" 1 km s~ 1 ) a 


120±?° 


80±H 


59 ±6 


25 ±3 


AI(v p ) (Jy beanr 1 ) 3 


6.5±1.0 


2.4±0.4 


4.3±0.4 


1.3±0.2 


ATi> p ) (K) a - b 


112±17 


41±7 


> 100 


> 30 


/e b 


0.6 


0.6 


< 1 


< 1 


T ox (K) 


120 


65 


120 


65 


TL(fp) 


3.9 


1.3 


1.9 


0.7 


Nh 2 (10 22 cm- 2 ) c 


6.4 


1.3 






Mout (M ) 
tout (yr) d 


1.3 


0.1 






1600/tan(i) 


500/tan(i) 







Notes: 

a The quantities at the maximum positions of the CO(3-2) flux. AV = J AI(v)&v / AI(v p ). 

b The beam filling factor (/b) is the ratio of the minor-axis size to the beam size projected on the minor axis and 
AT L (/, P ) = ^V 2 A/K)/ B 1 . 

°JVh 2 = 8 x 10 17 T cx exp [l - exp(^^)] ' / r L dv « 5 x 10 16 cm- 2 T 2 x exp (^f) t l (u p )AV for CO(3-2) 

assuming an abundance of [iVco] / [Nh 2 ] — 10~ 4 . 
i is the angle between the major axis of outflow and line of sight. 



which corresponds to r,; < 0.03 pc from the resu lt of 
3.6 cm continuum observation bv iMehringerJ (|1994f h We 
derived the Stromgren radius of 0.008 pc using the mean 
physical parameters inferred above. The sound speed of 
the ionized medium is C\\ — v/ jkT e /m p , where m p is 
the mass of the proton and k is the Boltzmann constant. 
Assuming the adiabatic index 7 = § and an average 
T e = 4900K, we found C n « 8.2km s _1 . Therefore, the 
dynamic age for the ultracompact Hn region W51e2-W 
is about 5000 yr, which appears to be a few times older 
than the outflow from W51e2-E (see t out in Table [T]). 

High-resolution observations of absorption spectra of 
the HCN(4-3) from W51e2-E and W51e2-W show that 
at present W51e2-E dominates the gas accretion from 
the molecular core W51e2 (see the HCN(4-3) absorption 
region, blue color in Figure [5]). If the accretion rate is a 
constant of 1.3 xlO" 3 M Q yr" 3 (|Shi et al.ll201fl . W51e2- 
E would have accreted a total of ~6 M Q over the 5000 yr 
period, <5% of the total mass (140 M Q ; Shi et al. 2010) 
of the sub-core. Thus, most of the mass accumulation for 
W51e2-E probably occurred before the star(s) formed in 
the ultracompact Hn region W51e2-W. After an 08 star 
or a cluster of B type stars formed in W51e2-W, the 
overwhelming radiation pressure and strong turbulence 
in W51e2-W halted the accretion in the sub-core region. 
Therefore, W51e2-E might speed up its formation of stars 
since W51e2-W is no longer competing for accretion. 

In summary, a bipolar molecular outflow in W51e2 is 
confirmed and the massive protostellar core W51e2-E is 
identified to be the origin of this powerful outflow. In 
comparison to the ultracompact Hn region W51e2-W, 
W51e2-E appears to be the dominant accretion source 




19 h 23 m 44.10 s 44.00 s 43.90 s 43.80 s 

Right Ascension (J2000) 



Fig. 2. — RGB composition image of H26a emission (red), 0.85 
mm continuum emission (green), and HCN(4-3) absorption (blue) 
observed with t he SMA at a resolution of 0."33 X 0."24 (P.A.=68°; 
IShi et aL1l2010r ) overlaid with the contour map of CO (3-2) outflow 
(red and blue) and absorption (grey) (Figure [l^). Crosses mark 
the positions of the sub-cores. 



at present, where the active formation of massive stars 
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takes place. 
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